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ABSTRACT 



> 

, We use hydrodynamical simulations of galaxy clusters and groups to study the effect 

■ of pre-heating on the entropy structure of the intra-cluster medium (ICM). Our simu- 
' lations account for non-gravitational heating of the gas either by imposing a minimum 

l/^ , entropy floor at redshift z/j = 3 in adiabatic simulations, or by considering feedback 

' by galactic winds powered by supernova (SN) energy in runs that include radiative 

cooling and star formation. In the adiabatic simulations we find that the entropy is 
" ^1 increased out to the external regions of the simulated halos as a consequence of the 

I ' transition from clumpy to smooth accretion induced by extra heating. This result is in 

O I line with the predictions of the semi-analytical model by Voit et al. However, the in- 

troduction of radiative cooling substantially reduces this entropy amplification effect. 
^ ' While we find that galactic winds of increasing strength are effective in regulating star 

formation, they have a negligible effect on the entropy profile of cluster-sized halos. 
Only in models where the action of the winds is complemented with diffuse heating 
corresponding to a pre-coUapse entropy do we find a sizable entropy amplification out 
, to the virial radius of the groups. Observational evidence for entropy amplification in 

■ the outskirts of galaxy clusters and groups therefore favours a scenario for feedback 

that distributes heating energy in a more diffuse way than predicted by the model for 
galactic winds from SN explosions explored here. 

Key words: Cosmology: numerical simulations - galaxies: clusters - hydrodynamics 
- X-ray: galaxies 



1 INTRODUCTION 

The thermodynamic properties of the intra-cluster medium 
(ICM) reflect the full history of the cosmic evolution of 
baryons. This evolution is primarily driven by gravity, which 
shapes the gas distribution on large scales and determines 
the hierarchical assembly of cosmic structures, but also by a 
complex interplay with the hydrodynamic processes of star 
formation and galaxy evolution (see Rosati et al. 2002; Voit 
2005 for recent reviews, and references therein). Due to the 
high density and temperature reached by gas in the envi- 
ronment of galaxy clusters and groups, X-ray observations 



provided so far the most powerful diagnostics on its thermal 
properties. 

The simplest model for the evolution of the ICM ac- 
counts only for the action of gravity (Kaiser 1986). This 
picture is said to be self-similar because it predicts that 
galaxy groups and clusters of different richness are essen- 
tially scaled versions of each other. Based on this realization, 
one can then derive unique scaling relations between X-ray 
observables. For instance, the X-ray luminosity, Lx , should 
scale with the gas temperature T as Lx oc T^; the "entropy" 
S — T jn'J'^ (rie is the electron number density) should scale 
linearly with temperature owing to the self-similarity of the 
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gas density profiles. Furtliermore, splierical gas accretion in 
a NFW halo (Navarro, Frenk & White 1997) predicts en- 
tropy to scale with the cluster-centric distance as S oc R^'^ 
(Tozzi & Norman 2001). 

While adiabatic hydrodynamical simulation of clusters 
have validated this simple picture (e.g., Navarro et al. 1995; 
Eke et al. 1998; Borgani et al. 2001), a number of obser- 
vations have now established that gravity cannot be the 
only player: the gas density profiles in the central regions 
of groups and poor clusters are observed to be shallower 
than in the self-similar model, and the relative entropy level 
is correspondingly higher than in rich clusters (e.g., Ponman 
et al. 1999; Finogucnov ct al. 2002; Ponman ot al. 2003). As a 
consequence, smaller systems are relatively undcrluminous, 
thus producing too steep an Lx-T relation, Lx oc T^^ (e.g., 
Markevitch 1998; Arnaud & Evrard 1999; Sanderson et al. 
2003). 

Those observational facts call for additional physics in 
the form of non-gravitational processes, capable of break- 
ing the self-similarity of the ICM properties. As a solution, 
a number of authors have proposed the presence of extra 
heating, possibly occurring before the cluster collapse (e.g. 
Evrard & Henry 1991; Kaiser 1991; Bower 1997; Cavaliere, 
Menci & Tozzi 1998; Balogh et al. 1999; Bialek, Evrard & 
Mohr 2001; Tozzi & Norman 2001; Borgani et al. 2002; Dos 
Santos & Dore 2002), or of radiative cooling (e.g., Bryan 
2000; Voit & Bryan 2001; Wu & Xue 2002), or, perhaps more 
likely, of a combination of these two effects (e.g., Voit ct al. 
2002; Muanwong et al. 2002; Tornatore et al. 2003; Borgani 
et al. 2004; Kay et al. 2003, 2004). Radiative cooling deter- 
mines the minimum entropy level of the X-ray emitting gas 
through a selective removal of gas with short cooling time. 
On the other hand, non-gravitational heating affects the 
amount of gas which can cool down to form stars, thereby 
preventing a cooling runaway. Independent of the details 
of the model and the nature of the assumed heating en- 
ergy, it has been common prejudice that the ICM entropy 
would only be affected in the central regions of clusters, 
while the self-similarity would be preserved in the halo out- 
skirts, where dynamics is still dominated by gravity. 

However, based on a semi- analytic model for gas accre- 
tion in clusters, Voit et al. (2003) have recently predicted 
that a smoothing of the gas density due to preheating in 
infalling sub-halos would boost the entropy production at 
the accretion shock of clusters. As a result, an excess of 
entropy with respect to the prediction of the self-similar 
model is generated in the cluster outskirts. This effect due to 
smooth accretion should be more important for poorer sys- 
tems, since their accretion is dominated by smaller clumps 
which are more affected by pre-heating. Therefore, entropy 
generation is expected to be more efficient in smaller sys- 
tems, thus boosting their entropy profiles more than those 
of more massive structures. 

Observational support for this model has been provided 
by Ponman et al. (2003). Using ROSAT/ASCA data for a 
large set of clusters and groups, they found that the entropy 
level in the outskirts is higher than predicted by self-similar 
scaling, in fact by a larger amount for poorer systems. Pon- 
man et al. (2003) also gave a possible interpretation in terms 
of smoothing of the accretion pattern, concentrating on the 
effect of smoothing on accreting filaments rather than on gas 
clumps. The lack of self-similar scaling for the overall nor- 



malisation of the entropy profiles has also been confirmed by 
Pratt & Arnaud (2004) and Piffaretti et al. (2005), based on 
XMM-Newton observations, and by Mushotzky et al. (2003) 
from Chandra data. Thanks to the possibility to perform 
spatially resolved spectroscopy, these authors were able to 
derive the entropy profiles from the actual temperature pro- 
files. Consistent with the results by Ponman et al. (2003), 
they concluded that the shape of the entropy profiles ap- 
pears to be quite similar for all systems, with an amplitude 
that follows the scaling S oc T^/"* more closely than the 
S ocT scaling expected from self-similarity. 

The question then arises as to whether the observed en- 
tropy structure of groups and clusters can be used to derive 
information on the nature of the feedback that affects the 
thermodynamical history of cosmic baryous. Whatever the 
astrophysical source of this feedback is, it must act in such 
a way as to regulate star formation and, at the same time, 
to reproduce the observed entropy profiles. 

Based on the model of star formation and supernova 
(SN) feedback by Springel & Hernquist (2003a, SH03 here- 
after), Borgani et al. (2004) performed a large cosmologi- 
cal hydrodynamical simulation with the aim of studying the 
X-ray properties of galaxy clusters and groups. Although 
the model provides a realistic description of the cosmic star 
formation history (Springel & Hernquist 2003b), clusters 
are found to have approximately self-similar entropy pro- 
files, which is not consistent with observations. This result 
demonstrates that generating the required level of entropy 
amplification at cluster outskirts is a challenging t£isk for 
hydrodynamical simulations even if they treat star forma- 
tion and feedback. A number of simple feedback schemes 
have been suggested that heat the gas surrounding galaxies 
in order to better reproduce the observed entropy scaling 
(e.g. Kay et al. 2004). While such experiments provide very 
useful guidelines, it is clear that a self-consistent numeri- 
cal implementation of a physically well motivated feedback 
model which can successfully satisfy a large body of obser- 
vatioual constraints is highly desirable. 

In this paper, we use hydrodynamical simulations to 
explore a suite of different heatiug schemes, both for radia- 
tive and non-radiative runs. We focus on four objects which 
cover the mass range from poor clusters to poor groups. The 
aim of our analysis is to verify whether an appreciable en- 
tropy amplification at large radii, as predicted by Voit et 
al. (2003), can be obtained with simple preheating schemes 
and/or with self-consistent models that include feedback 
from SN-driven galactic winds. 

The main questions that we intend to address with our 
analysis are the following: (a) Is the entropy amplification 
by smoothed accretion confirmed by hydrodynamical simu- 
lations when a simple scenario for pre-hcating with an en- 
tropy floor is invoked? (b) Does this effect persist when sim- 
ulations include a realistic treatment for star formation and 
SN feedback? (c) What are the observed cutropy profiles 
telling us about the nature of the feedback that affects the 
evolution of the ICM and IGM? 

The outline of this paper is as follows. In Section 2, we 
describe the simulations that we carried out and the differ- 
ent schemes to provide non gravitational heating. Section 3 
is devoted to the presentation and discussion of the results, 
both for the non-radiative and for the radiative runs. In Sec- 
tion 4, we summarise our results and give our conclusions. 
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2 THE SIMULATIONS 

Our simulations have been carried out with GADGET-2, 
a new version of the parallel Tree-SPH code GADGET 
(Springel, Yoshida & White 2001). It uses an entropy- 
conserving formulation of SPH (Springel & Hernquist 2002), 
and includes radiative cooling/heating by a uniform evolving 
UV background, a sub-resolution description of star forma- 
tion from a multiphase interstellar medium, and the effect 
of galactic winds powered by SN energy feedback (SH03). 

The simulated groups and clusters are the same as in 
Tornatore et al. (2003) and Finoguenov et al. (2003), and 
we refer to these papers for more details on their char- 
acteristics. Our set of simulated objects includes a mod- 
erately rich cluster with Mvir ^ 2.6 x lO^'^/i^^M©, and 
three groups-sized objects having virial mass in the range 
1.6 X 10"< Mvir/(/i"^M0)< 4.2 X 10". The correspond- 
ing halos have been selected from a DM-only simulation of 
a ACDM model with Qm = 0.3, h = 0.7, as = 0.8 and 
/bar = 0.13 for the baryon fraction, within a box of size 
70/i~^Mpc. As such, they encompass the mass range char- 
acteristic of moderately rich clusters to groups, where the 
effect of non-gravitational heating is expected to be impor- 
tant. 

Using a re-simulation technique, the mass and force res- 
olution has been increased in the Lagrangian regions of the 
parent box which correspond by a = to volumes encom- 
passing several virial radii around each of the selected halos. 
The resolution is progressively degraded in regions farther 
away, allowing us to save computing time while still provid- 
ing an accurate representation of the large-scale tidal field. 
In the high resolution region of the cluster simulation, gas 
particles have mass mgas — 2.2 x W^h^^Mr.), and wc chose 
epi = 5h~^ kpc for the Plummer-equivalent softening scale, 
kept fixed in physical units out to « = 2, and then fixed in 
co-moving units at earlier times. In order to have a com- 
parable number of particles within the virial regions of the 
different simulated structures, we increased the mass reso- 
lution by a factor of eight for the runs of the groups and, 
correspondingly, we decreased the force softening by a factor 
of two. The main characteristics of the simulations are listed 
in Table 1. 

We have simulated each object both with radiative and 
non-radiative physics. Our non-radiative runs include the 
case of purely gravitational heating (GH) and two pre- 
heating models, which are implemented by imposing a min- 
imum gas entropy of Sa = 25 and 100 keV cm^ at redshift 
Zh = 3, respectively. Entropy of heated particles is raised 
by incresing their internal energy (i.e., temperature), while 
leaving their density unchanged. For the radiative runs, we 
have performed several simulations by changing either the 
parameters controlling the galactic winds or the level of the 
pre-heating entropy floor. In the model of SH03, the ve- 
locity of the galactic winds, v^, scales with the fraction rj 
of the SN-II feedback energy that contributes to the winds, 
as Vw oc rj^^'^ (see eq.[28] in SH03). The total energy pro- 
vided by SN-II is computed by assuming that they are due 
to exploding stars with meiss > 8 Mq , each SN releasing 
10'^'^ ergs, and that their abundance is given by a Salpctcr 
(1955) initial mass function (IMF). In our simulations, we 
consider the cases 77 = 0.5, 1 and 3, yielding Vw ~ 340, 480 
and 840 km s~^ (runs Wl, W2 and W3), respectively. While 



Table 1 . Physical characteristics and numerical parameters of the 
simulated halos. Column 2: total mass within the virial radius at 
z = (lOi^/i~^M0); Column 3: virial radius ( h~^Mpcy, Column 
4: mass-weighted temperature within Rvir (keV); Column 5: mass 

of the gas particles in the simulations (10*/i~"'^Mq): Column 6: 
Plummer-equivalent gravitational softening at z = (fe^-'^kpc). 
The numbers reported in Columns 2-4 refer to the non-radiative 
run with gravitational heating (GH) only, but their values do not 
significantly change for the other runs. 



Run 




-Rvir 


^ vir 




e 


Cluster 


26.4 


1.31 


2.04 


2.17 


5.0 


Group- 1 


4.17 


0.71 


0.52 


0.27 


2.5 


Group-2 


1.75 


0.53 


0.39 


0.27 


2.5 


Group- 3 


1.64 


0.52 


0.28 


0.27 


2.5 



assuming r/ > 1 is clearly unrealistic in this picture, it can 
be phenomenologically interpreted as being due to an extra 
energy source, e.g. coming from a top-heavier IMF or from 
an AGN component. 

In their model for galactic winds, SII03 treated SPH 
particles that become part of the wind as temporarily de- 
coupled from hydrodynamical interactions, in order to allow 
the wind particles to leave the dense interstellar medium 
without disrupting it. This decoupling is regulated by two 
parameters. The first parameter, pdcc, defines the minimum 
density the wind particles can reach before being coupled 
again. If pth is the threshold gas density for the onset of star 
formation (pth ^ 2.8 x 10"^^g cm"^ in the SH03 model), 
then it should be pdcc < Pth for the winds to leave the star- 
forming region, and the typical setting by SH03 has been 
Pdec = 0.1 — 1.0 Pth. The second parameter, Zdec, provides a 
maximum time via tdec = Idec/ve a wind particle may travel 
freely before becoming hydrodynamically coupled again. If 
this time has elapsed, the particle is coupled again, even if 
it has not yet reached pdcc. The /dec parameter is introduced 
only to ensure that the decoupling is stopped quickly also in 
cases where the wind particle cannot escape the dense ISM 
due to gravity alone. We provide the values of pdcc and /doc 
used for the radiative runs in Table 2. The Wl, W2 and W3 
models all have the same values for these parameters, but 
they differ in their wind velocity. In an attempt to produce 
more wodely distributed by winds, we have also run an addi- 
tional model W4 which has the same wind speed as W3, but 
where we allowed the winds to reach much lower densities 
and larger distances while being decoupled. 

In summary, we have performed 9 runs for each of the 
four simulated structures, corresponding to different ways of 
changing the gas physics. The characteristics of these runs 
are given in Table 2. We also provide in this table the mean 
value of non-gravitational specific energy assigned to the 
gas that falls at the end of the run within the virial radius 
of the cluster (comparable amounts are found for the three 
groups). For the simulations with pre-heating, wc have de- 
termined this energy by tracking back to 2: = 3 all the gas 
particles within R^h at « = and summing up the energy 
contributions required to increase their entropy to the floor 
value where needed. For the runs including SN feedback, 
we measure the total stellar mass within R^ir aX z = and 
we compute the total number of SN-II expected for the as- 
sumed Salpeter IMF. The resulting energy is then multiplied 



© 2003 RAS, MNRAS 000, 000-000 



4 S. Borgani et al. 




Figure 1. Entropy profiles of the non— radiative simulations. In the left and in the right panels, we compare the results for the 
gravitational-heating (GH) runs (solid curves) to those for the pre-heated runs (short-dashed curves) with entropy floor = 25 
keV cm^ (S25) and 5fl = 100 keV cm^ (SlOO), respectively (see Table 2). Upper and lower panels correspond to the simulations of 
the Cluster and of the Group-3. In each panel, we also show with a long-dashed curve the entropy profile obtained by adding the 
entropy— floor value, Sfi, to the GH case. The difference between the short-dashed and the long-dashed curves provides the net entropy 
amplification. The two vertical dotted lines mark the positions of -R200 and -R50O1 defined as the radii encompassing an average density 
200pcr and SOOpcr, respectively (here pcr is the critical cosmic density). 



Table 2. Description of the different runs. Column 1: run name; 
Column 2: level of the entropy floor at 2^ = 3 (keV cm^); Col- 
umn 3: wind speed (km s~^); Column 4: specific heating en- 
ergy assigned to the gas particles that fall within i?vir at 2 = 
(keV/particle); Column 5: the limiting density for wind decou- 
pling (in units of the threshold density, pt^ , for star formation (see 
text); Column 6: the maximum length that winds travel while be- 
ing decoupled (units of /i~^kpc); Column 7: fraction of baryons 
in stars within iJvir- The values in Columns 4 and 7 refer to the 
simulations of the Cluster. 



Run 


Sa 




Eh 


Pdcc 


^dcc 




Non-rad. runs 














GH 














S25 


25 




0.5 








SlOO 


100 




2.2 








Rad. runs 














Wl 




341 


0.3 


0.5 


10 


0.20 


W2 




484 


0.5 


0.5 


10 


0.17 


W3 




837 


1.2 


0.5 


10 


0.14 


W4 




837 


1.1 


0.01 


50 


0.12 


W1+S25 


25 


341 


0.2-|-0.4t 


0.5 


10 


0.13 


WH-SIOO 


100 


341 


0.2-fl.8t 


0.5 


10 


0.13 



t We report separately the energy budget associated to galactic 
winds (see text) and the pre— heating energy required at 2 = 3 to 
establish the entropy floor. 



by the efficiency parameter i] to compute the total energy 
associated with wind feedback. 

Looking at Table 2, we note that the heating energy 
required to create the entropy floors is comparable to those 
of the winds. However, the way in which these two feedback 



schemes affect the thermodynamics of the diffuse baryons 
is intrinsically very different. The first heating mechanism 
provides a diffuse energy input that can in principle reach 
all the gas particles in the simulations, and is provided in 
an impulsive way (i.e. all the energy is released at the heat- 
ing redshift Zh = 3). The second mechanism, instead, re- 
leases energy gradually in time, since it follows the pattern 
of star formation, and it is narrowly concentrated around 
star-forming regions. As we will discuss in the following, 
these differences have an important impact on the smooth- 
ing of the accretion pattern and, therefore, on the entropy 
generation by accretion shocks. 



3 RESULTS 

3.1 Non— radiative runs 

In Figure 1, we compare the effect of imposing the two dif- 
ferent entropy fioors on the entropy profiles of the cluster 
simulation and of 'Group-3'. Here and in the following we 
plot profiles to an innermost radius which contains 100 SPH 
particles. This radius has been shown to be the smallest one 
where numerically converged results for the X-ray luminos- 
ity can be obtained (Borgani et al. 2002). 

As expected, the effect of pre-heating is that of increas- 
ing the level of the ICM entropy by a larger amount for the 
less massive system, an effect that increases in strength for 
a higher entropy fioor. The runs with gravitational heating 
display a number of wiggles which mark the positions of 
merging sub-halos carrying low-entropy gas. These wiggles 
are erased once the gas is pre-heated. This is consistent with 
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the expectation that pre-heating destroys the gas content of 
halos whose virial temperature is lower than the specific pre- 
heating energy. In the cluster run, the entropy increase due 
to pre-heating is marginal at the virial radius, thus implying 
that the overall ICM thermal energy content within Jivir is 
mainly determined by the action of gravity. 

In Figure 2, we provide further information about the ef- 
fect of pre-heating on the entropy structure of the ICM. Here 
we show the entropy-density phase diagram of gas particles. 
The two separate concentrations of points mark the gas par- 
ticles falling within 0.1 i?vir and those having cluster-centric 
distances in the range (0.9—1.1) i?vir, the former being those 
having larger densities. Following these two populations of 
gas particles allows us to understand the different effects 
that non-gravitational physics has in the innermost and in 
the outskirt regions of the ICM. In the outer cluster regions 
in the GH runs (left panels), the change of slope in the S-pg 
relation marks the transition from diffuse (shocked) accre- 
tion to clumpy accretion (see also Sect. 3.2 and Figure 5 
below). The clumpy accretion is associated with merging 
subgroups and filaments, which are characterised by low en- 
tropy and high density, and are erased by pre-heating with 
an entropy floor of 100 keV cm^ (SlOO runs, right panels). 
This demonstrates that the effects of pre-heating are: (a) 
to cancel any signature of clumpy accretion in the outskirts; 
(b) to suppress the gas density in the central halo regions, 
and (c) to amplify the entropy of both the diffuse phase in 
the outer regions and of the gas in the central regions, by 
an amount that increases for lower-mass systems. 

In their semi-analytic description of entropy generation 
from diffuse accretion, Voit et al. (2003) have shown that 
the post-shock gas entropy is higher than it would be with- 
out preheating, at a given density, by an amount between 
0.84 S and S, where S is the added entropy. Therefore, when 
comparing the profiles of preheated simulations to the cor- 
responding GH profiles where the preheating entropy level 
is simply added to it, any excess entropy should then bo 
attributed to the amplification effect due to accretion of 
smoothed gas instead of clumps. We expect this effect to be- 
come more prominent for stronger pre-heating, which makes 
the accretion more diffuse, and for smaller groups, where the 
accretion is dominated by comparatively smaller gas clumps. 
In Fig. 1, the long-dashed curves show the profiles obtained 
by adding the floor value, Sb, to the GH profile. Consistent 
with the expectation from the model by Voit et al. (2003), 
we find clear evidence for the entropy amplification effect. 
It is more apparent for the higher entropy floor and for the 
smaller halo of the Group-3. For instance, at i? = 0.1i?vir, 
we obtain for the cluster simulations an amplification of 14 
per cent and 84 per cent for the S25 and SlOO runs, respec- 
tively. These numbers increase to about 230 per cent and 
300 per cent, respectively, for the 'Group-3' runs. 

We note that our pre heating scheme does not increase 
the entropy of each gas particle by a fixed extra amount. 
Instead, it creates an entropy fioor, so that the extra entropy 
assigned to each particle can be zero if the particle is already 
at high entropy, or just the value required to attain the floor 
value. On the other hand, the model by Voit et al. (2003) is a 
prediction for the entropy amplification when pre-heating is 
implemented by adding a constant amount of entropy to all 
gas particles. In this case, we would have obtained a stronger 



pre-heating and, therefore, an even stronger amplification 
effect. 

The signature of differential entropy amplification ex- 
tending out to outer regions of groups and clusters is in line 
with observational evidences based both on ROSAT/ASCA 
(e.g. Ponman et al. 2003) and XMM-Newton data (e.g. 
Pratt & Arnaud 2004; Piffaretti et al. 2005). However, de- 
spite the trend of progressively higher relative entropy lev- 
els in groups, the observational data also indicates that the 
shape of entropy profiles is almost independent of mass, with 
no evidence for the presence of large isentropic cores. This is 
quite different from the behaviour seen in our non-radiative 
simulations, where an amplification of entropy at the clus- 
ter outskirts is obtained at the price of creating much flatter 
entropy profiles in groups and larger isentropic cores than 
in clusters (see Fig.l). 

3.2 Simulations with radiative cooling 

The results for the entropy profiles change substantially for 
the radiative runs. In Figure 3, we compare the profiles for 
the simulations including cooling and star formation, ob- 
tained by either changing the parameters of the winds (left 
panels) or by adding a pre-heating entropy fioor to the feed- 
back of the Wl model. As for the "Cluster" , increasing the 
wind speed leaves the entropy profiles completely unaffected 
in the halo outskirts. As for the central regions, increasing 
the wind speed from the Wl to the W2 model slightly de- 
creases the entropy level. In fact, the stronger feedback cre- 
ates a population of lower entropy particles which are kept 
in the hot phase by the larger amount of energy feeback, 
instead of cooling down. Further increasing the wind speed 
(W2 and W3 runs) or changing the parameters defining the 
wind decoupling do not further change the entropy profile 
for the Cluster. 

Although the effect of stronger winds is more appar- 
ent for the 'Group-3' run, it is nevertheless still marginal 
for Rf^ 0.5ifvir. However, increasing the wind speed has a 
significant effect on the entropy profile in the inner regions 
of 'Group-3': the stronger feedback creates a nearly isen- 
tropic regime for 0.07 ^ R/Rvii^ 0.5, while the profile steep- 
ens again in the innermost regions, where cooling starts to 
dominate. Quite intriguingly, a similar behaviour has been 
recently found by Mahdavi et al. (2005), based on an analy- 
sis of galaxy groups observed with XMM-Newton. They find 
entropy profiles with a broken power-law shape S oc R", 
with an inner and outer slope of ~ 0.9 and ~ 0.4, respec- 
tively, and a transition that takes place at about 0.1J?5oo- 
We also note that the presence of a strong wind reduces the 
resulting fraction of stars from about 20 per cent to about 13 
per cent (see Table 2). On one hand this demonstrates that 
strong winds do increase entropy and, therefore, the cool- 
ing time of the gas surrounding star-forming regions (e.g., 
Springel & Hernquist 2003b; Borgani et al. 2004, in prepa- 
ration). On the other hand, the results shown in the left 
panels of Fig. 3 also indicate that this effect is not causing 
a transition from clumpy to diffuse accretion, as would be 
required to generate an entropy amplification effect. 

Allowing the wind particles to be decoupled until they 
reach lower densities, and hence allowing them to travel to 
a larger distance from the star-forming regions before in- 
teracting with the diffuse gas (model W4), only produces a 
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modest change. Since the mechanical wind energy is ther- 
malized at lower gas densities in this model, radiative losses 
are reduced as a consequence of longer cooling times. How- 
ever, the resulting additional suppression of star fraction 
(see Table 2) and increase of the entropy in the 'Group-3' 
simulation are too small to change our conclusions. 

These results on the effect of winds are also supported 
by the phase-diagrams shown in Figure 4, by the relative 
distribution of low-entropy and high-entropy particles in 
the halo outskirts, shown in Figure 5, and by the maps of 
gas density of the 'Cluster' runs, shown in Figure 6. The 
latter demonstrates that both in the proto-cluster region at 
z = 2 and in the cluster at z = 0, the stronger winds wash 
out only the smallest halos and make the larger ones slightly 
pufBer, while preserving the general structure of the cosmic 
web surrounding the Lagrangian cluster region. 

By comparing the results of the Wl and W3 models in 
Figs. 4 and 5 (left and central panels, respectively), we see 
that a higher wind speed has also a negligible effect on the 
Cluster, despite being able to erase accreting clumps in the 
outskirts and to push the gas in the central regions of the 



'Group-3' simulation to higher entropy. In all cases, we find 
that the stronger winds of the W3 and W4 models are quite 
ineffective in increasing the entropy of the diffusely accreting 
gas. 

A larger effect on the ICM entropy structure is induced 
by adding a pre-heating entropy floor to the Wl feedback 
(right panels of Fig. 3). In this case, a sizable amplification 
is obtained for both the 'Group-3' and the 'Cluster' simula- 
tions, at least when the higher entropy fioor, Sn — 100 keV 
cm^, is used. As shown in Fig. 6, this entropy amplification 
is associated with a much smoother gas density distribution, 
both at 2 = 2 and at z = 0. In this case, the filamentary 
structure of the gas distribution is completely erased, while 
only the largest halos are able to retain part of their gas 
content. This visually demonstrates that a transition from 
clumpy to smooth accretion during the process of cluster 
formation requires a more broadly distributed form of feed- 
back. 

The effect on the S-pg phase diagram is shown in the 
right panels of Fig. 4: Adding the entropy ffoor has a stronger 
impact in smoothing the pattern of gas accretion than in- 
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Figure 5. The distribution of gas particles in the outskirts (0.9 < R/Rvir < 1.1) of the Cluster (upper panel) and Group-3 (lower panels). 
Left, central and right panels are for the Wl, W3 and Wl-l-SlOO models, respectively. Light and heavy points are for the gas particles 

whose entropy is above and below the value of 400 keV cm^ for the Cluster and of 60 kcV cm^ for the Group-3. This plot demonstrates 
that such entropy values approximately mark the transition between clumpy and diffuse gas for the two simulated structures. 



creasing the wind speed. Also, Fig. 5 shows that almost no 
signature of clumpy accretion is left in the 'Group-3' run. 
Comparing these results to those of the non-radiative runs 
shows that radiative cooling significantly reduces the en- 
tropy amplification. The main reason for this is that cooling 
increases gas density within accreting clumps and filaments, 
thus partially inhibiting the transition from clumpy to dif- 
fuse accretion. This effect can be appreciated by comparing 
the left panels of Fig. 4 with the left panels of Fig. 2. While 
the presence of cooling does not significantly change the en- 
tropy level of the accreting difi^use gas, it nevertheless allows 
residual low -entropy clumps in the outskirts to survive the 
pre-heating, thus lowering the overall amplitude of the en- 
tropy profiles. 

An interesting characterisation of the effect of either 
increasing the wind speed or adding an entropy floor is ob- 
tained by comparing the star fraction with the correspond- 
ing specific extra energy involved by non-gravitational heat- 
ing (see Table 2). Increasing the wind speed in the 'Cluster' 
runs has only a marginal effect on the entropy profiles, but 
it is quite effective in regulating cooling: the W3 model pro- 
duces ~ 30 per cent fewer stars than the Wl model, bringing 
it into better agreement with observational data (e.g., Lin, 



Mohr & Stanford 2003) which generally favour low stellar 

densities. Quite interestingly, strong winds arc at least as 
efficient as the SlOO model in preventing overcooling, while 
requiring a lower amount of extra energy per gas particle. 
This highlights the different ways in which winds and an en- 
tropy floor affect the thermodynamics of the diffuse baryons. 
The former acts locally to increase the cooling time of the 
gas surrounding the star forming regions but it is inefficient 
in smoothing out the gas content of merging sub-halos. The 
latter provides a diffuse impulsive heating and provides for 
a much more efficient transition from smooth to diffuse ac- 
cretion. 

A further question we investigated is whether the dif- 
ferent feedback and pre-heating schemes are able to violate 
the self-similarity of the entropy profiles as strongly as it is 
observed. For this purpose, we plot in Figure 7 the profiles of 
reduced entropy, S'(r)/T'mw(f), for the four simulated struc- 
tures in different radiative runs. Since in this analysis we do 
not aim for a detailed comparison with observations, we here 
prefer to normalise the entropy to the mass weighted tem- 
perature, Tmw, instead of using the emission-weighted tem- 
perature. The latter has been shown to provide an inaccurate 
measure for the spectroscopic temperature of the ICM that 
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Figure 6. Maps of the gas density for the radiative runs of the Cluster. Upper, central and lower panels are for the Wl, W3 and 
Wl+SlOO runs, respectively. Left and right panels are for the outputs at 2 = 2 and 2 = 0, respectively. At 2 = the size of the box is 
11.7/i~^Mpc, while at 2 = 2 it corresponds to 17.5?t~^Mpc comoving. The small white knots mark the "galaxies", i.e. the places where 
high density gas is undergoing cooling and star formation. 
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Figure 7. The profile of reduced entropy for the radiative runs. The left panels and the upper right panel show the results for the runs 

with wind velocity v-\jv = 341kms~^ (Wl runs) and 837kms~^ (W4 runs), and for the Wl runs that also include an entropy floor of 100 

keV cm^ (Wl+SlOO), respectively. The profiles are scaled with the mass— weighted temperature, Tmw. The lower right panel shows the 

2/3 

results of the Wl+SlOO runs, but with entropy profiles normalised according to the observed temperature scaling, . In all panels, 
the straight dotted line marks the slope S oc RP '^^ which is the best— fit to the observed entropy profiles of Pratt & Arnaud (2004) and 
Piffaretti et al. (2005). The vertical dotted lines are for r2oo and rsoo of the Cluster. 



is observationally inferred from X-ray spectra (e.g., Math- 
iesen & Evrard 2001; Mazzotta et al. 2004; Vikhlinin 2005). 
Although cooling, star formation and SN feedback are likely 
to introduce characteristic scales in the ICM thermodynam- 
ics, all the Wl runs have remarkably self-similar entropy 
profiles, over the whole range of resolved scales. However, 
for the W4 runs, the self-similarity is clearly broken in the 
central regions, while it is largely preserved at Rfi O.Siivir- 
Over the range of scales where self-similarity is preserved, 
the slope of the entropy profiles is always consistent with 
the S oc behaviour found by Pratt & Arnaud (2004) 

and Piffaretti et al. (2005) from their analysis of groups and 
clusters observed with XMM-Newton. 

As expected, only the SlOO entropy fioor generates an 
appreciable breaking of self-similarity in the outer regions 
of the simulated structures. In the lower right panel of Fig- 
ure 7 we show the entropy profiles rescaled according to 
Tmvf , which is consistent with the scaling suggested by ob- 
servational data. Quite ironically, while this time we are able 
to reproduce the scaling with temperature of the normali- 
sation of the profiles, we lose the agreement with the radial 
slope, which is significantly shallower than the observed one. 
We note here that an accurate derivation of the ICM en- 
tropy profiles from observations requires properly resolving 
both the gas density profiles and the temperature profiles. 
While different authors generally agree on the radial run 
of gas density, the temperature profiles are still a matter 
of vigorous debate. For instance, both ASCA (e.g., Marke- 
vitch et al. 1998; Finoguenov et al. 2001) and Beppo-SAX 
(e.g., De Grandi & Molendi 2002) consistently find decreas- 
ing profiles, at least for 0.2i?iso (cf. also White 2000). 
However, Pratt & Arnaud (2003, 2004) find profiles consis- 
tent with being isothermal, at least outside the cooling re- 



gions. More recently, Piffaretti et al. (2005) and Vikhlinin et 
al. (2004) find decreasing outer profiles from XMM-Newton 
and Chandra data, respectively. While we do not intend to 
enter this debate here, we note that shallower entropy pro- 
files should be obtained from data when the corresponding 
temperature profiles have a negative outer gradient. 



4 DISCUSSION AND CONCLUSIONS 

In this paper, we have analysed hydrodynamical simula- 
tions of the formation of one galaxy cluster and three galaxy 
groups, with the main aim of studying the ability of different 
models of non-gravitational heating to cause an amplifica- 
tion of entropy generation by accretion shocks as a conse- 
quence of a transition from clumpy to smooth accretion (e.g. 
Voit & Ponman 2003). For this purpose, we have carried out 
non-radiative simulations as well as simulations with cooling 
and star formation, and combined those with two different 
models for non-gravitational heating. The extra heating has 
been supplied in two different forms: (a)hy imposing an en- 
tropy floor at = 3; (b) hy including the effect of galactic 
winds (Springel & Hernquist 2003a) of different strengths in 
the runs with star formation. 

The main results of our analysis give answers to the 
three questions posed in the Introduction, and can be sum- 
marised as follows. 

(a) Smoothing the accretion pattern by pre-heating in sim- 
ulations with non-radiative physics does indeed amplify the 
entropy generation out to the radius where accretion shocks 
are taking place. As predicted by semi-analytic models (Voit 
et al. 2003), this amplification is more pronounced for lower- 
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mass systems, because they accrete from smaller sub-halos 
whose gas content is more easily smoothed by extra heating. 

(b) Radiative cooling reduces this amplification efiiect by a 
significant amount. Cooling has the effect of increasing the 
dumpiness of the accretion pattern, and thus acts against 
the smoothing induced by extra heating. We find that only 
our pre-heating model with an entropy floor Sr = 100 keV 
cm^ generates a sizable entropy amplification out to the ex- 
ternal regions, 0.5-Rvir, of group-sized halos. 

(c) The heating from galactic winds is efficient in regulat- 
ing star formation and in providing an increase of the ICM 
entropy in the central regions of groups. However, the winds 
have a negligible effect in smoothing the accretion pattern 
and thus do not trigger an appreciable entropy amplification 
effect. This result, which holds also for the strongest winds 
considered in our analysis, indicates that the role played by 
these winds tends to be fairly well localised around the star 
forming regions and therefore hardly affects the gas dynam- 
ics over the whole interior of dark matter halos. 

Our results show that the temperature-entropy scaling 
provides powerful constraints on the physical nature of the 

energy feedback that affects the thermodynamic history of 
the diffuse baryons. Smoothing the gas accretion pattern 
within halos of forming galaxy clusters recjuircs a rather 
non-local feedback, whose action should not only be that 
of preventing gas overcooling. The difficulties faced by our 
model for galactic winds suggest that we may miss a proper 
description of physics which is able to distribute the SN en- 
ergy in the diffuse medium more efficiently (e.g. Kay et al. 
2003, 2004, where cold galactic gas was impulsively heated to 
high temperature) . Alternatively, this may be taken as an in- 
dication that other astrophysical sources of energy feedback 
are required. In this context, the most obvious candidate are 
probably AGN, whose effects should be taken into account 
self-consistently in future cosmological hydrodynamical sim- 
ulations. 

In this paper, we have not yet performed an in-depth 
comparison with observational data because we preferred to 
first provide a general interpretative framework for X-ray 
observations. However, a detailed comparison with observa- 
tions in future work is very promising for extracting yet more 
information from the ICM thermodynamics. Particularly the 
incrcEising amount of high quality X-ray observations at the 
group scale (e.g. Finoguenov et al. 2004) together with the 
ever more sophisticated simulation models give a real hope 
that the intra-cluster baryons will remain an extremely use- 
ful tracer and will eventually allow us to understand the 
nature of feedback in groups and clusters. 
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